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ABSTRACT 

We present a detailed multi-epoch analysis of 31 potential southern hemisphere ra- 
dio calibrators that were originally observed as part of a program to maintain the 
International Celestial Reference Frame (ICRF). At radio wavelengths, the primary 
calibrators are Active Galactic Nuclei (AGN), powerful radio emitters which exist at 
the centre of most galaxies. These are known to vary at all wavelengths at which 
they have been observed. By determining the amount of radio source structure and 
variability of these AGN, we determine their suitability as phase calibrators for long 
baseline radio interferometry at 2.3 GHz. For this purpose, we have used a set of com- 
plementary metrics to classify these 31 southern sources into five categories pertaining 
to their suitability as VLBI calibrators. We find that all of the sources in our sample 
would be good interferometric calibrators and almost ninety per cent would be very 
good calibrators. 

Key words: quasar: general - galaxies: jets - radio continuum: galaxies - techniques: 
interferometric 



1 INTRODUCTION 

High angular resolution observations of weak radio sources 
(where self calibration is not possible) require calibrator 
sources for correction of systematic effects and effects of 
the atmosphere on the measured visibilities. Atmospheric 
fluctuations cause perturbations in visibility phase which, if 
not corrected, seriously limit both the sensitivity and im- 
age quality of an interferometric array. Phase calibrators 
are also required for astrometric observations. An ideal cal- 
ibrator would look the same on all observing baselines. It 
should be bright, unresolved, or at least compact and should 
not vary. A calibrator source should also be separated from 
the target sources by as small an angle as possible in or- 
der to look along the same line of sight through the atmo- 
sphere. Therefore, it is desirable to have calibrator sources 
evenly distributed across the whole sky. In practice, at radio 
wavelengths, calibrators are mostly Active Galactic Nuclei 
(AGN), whose fundamental source of power is believed to 
be the accretion of matter onto a super-massive black hole 
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|eg. Reesll997l ') . They are known to vary at every wavelength 
at which they have been studied. AGN are very compact and 
isotropically distributed around the sky. They are also very 
distant objects and therefore, generally, have no discernible 
proper motions on the sky. It is these qualities that make 
them suitable as calibrators at radio wavelengths. 

Due to the limited number of radio telescopes, few sur- 
veys for calibra tors have been carried out in the sout h- 
ern h emisphere (|Oiha et al.ll2004al . [2OO5I : iFev et al.ll2004al lbl, 
2006). Hence, there are fewer known calibrators in the south. 
A major expansion of radio astronomy observing capability 
is underway in the southern hemisphere. Two SKA (Square 
Kilometre Array) precur sors, the Sou t h Afr ican MeerKAT 
(Karoo Array Telescope, iBooth et all d2009t)') a nd AS KAP 
(Australian SKA Pathfinder. Ijohnston et al.1 l|2008h ) are 
presently under construction, leading to the SKA itself. It is 
clear that interferometry and Very Long Baseline Interfer- 
ometry (VLBI) in the southern hemisphere needs a dense 
network of calibration sources at high resolutions and a 
range of frequencies. MeerKAT is the South African SKA 
demonstrator telescope. When completed, it will be the 
largest radio telescope array in the southern hemisphere. 
It will operate in two frequency ranges, 0.58-2.5 GHz and 
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8-14.5 GHZ (|Booth et al.ll2009f) . MeerKAT will participate 
in VLBI observations with the European VLBI Network and 
other VLBI arrays. The SKA will operate in the f requency 
range 0.7 to 25 GHz jCarilli fc Rawlingsl |2004| ; ISchilizzil 
120071 ) and is expected to have much longer baselines for 
which a southern VLBI calibrator list will be essential. 

The United States Naval Observatory, in collabora- 
tion with NASA, the Laboratoire dAstrophysique de Bor- 
deaux (LAB) and the National Radio Astronomy Observa- 
tory (NRAOJ) has, since 1994 been observing AGN, once 
every two months, through the Research and Development 
VLBI (RDV) experiments. This has led to a wealth of 
data of sources observed simultaneously at 2.3 GHz and 
8.4 GHz at regular epochs. Images of these AGN form the 
USNO's Radio Reference Frame Image Database (RRFID) 
and LAB's Bordeaux VLBI Image Database (BVIDp 
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ICollioud fc Charlotl 120091. The RRFID is a database of 
about 6700 images of over 700 AGN sources compiled from 
geodetic and astrometric VLBI experiments. The BVID con- 
tains over 1800 images of 824 radio sources. It is the goal 
of the ongoing RDV program to image the radio reference 
frame sources on a regular b asis, to moni t or the m for vari- 
ability and structural change. IPiner et al.l l|2007l ) discuss jet 
kinematics in a subset of the RRFID sources, concentrat- 
ing on the 8.4 GHz observations. Their survey is made up of 
RRFID sources that have been observed at 3 or more epochs 
from July 1994 to December 1998. 

In this paper, we focus on 2.3 GHz observations of a 
sample of 31 potential southern VLBI calibrators (with de- 
clinations between and —60°) from the kinematic survey of 
IPiner et al.l l|2007i ). The choice of 2.3 GHz was made because 
both MeerKAT and ASKAP will be operating at frequen- 
cies close to this. We seek to characterise this selection of 
southern radio sources, and determine their suitability as 
calibrators for southern VLBI experiments, especially those 
using MeerKAT and the SKA when completed. 

In the following sections, we will describe the observa- 
tions and data reduction process. We will then describe the 
imaging and modelfitting process and go on to analyse the 
data and discuss the results. 



2 OBSERVATIONS AND DATA REDUCTION 

Research and Development VLBI (RDV) experiments are 
carried out using the 10 antennas of the NRAO's VLBA 
and up to 10 other antennas across the globe, including 
Hartebeesthoek (South Africa) when available. The use of 
the global array greatly improves the u-v coverage. Harte- 
beesthoek greatly improves the u-v coverage for sources in 
the south, which make up the sample in this paper. 

The 31 southern sources included in our sample are 
listed in Table [T] along with their optical properties. The 
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Figure 1. Sky distribution of the southern sources plotted on an 
Aitoff equal-area projection of the celestial sphere. The dotted 
line represents the Galactic plane while the dashed line is the 
ecliptic. 

source sky distribution is shown in Figure [1] The RDV ob- 
serving epochs are from July 1994 to January 2008, giving 
a total of 32 epochs and an average of 20 epochs per source. 
About 100 sources were observed in each 24 hr observing 
run, with an average on-source time of 15 minutes. The on- 
source time is not continuous, but divided into scans of be- 
tween 1 and several minutes long, spaced in time to give 
optimal u-v coverage. 

Observations were made in a dual frequency bandwidth 
synthesis mode to facilitate delay measurements for astrom- 
etry. Observations in this mode also allow imaging at both 
frequency bands. Eight individual frequency bands (IFs) 
were recorded simultaneously, each 8 MHz wide, with 4 at 
2.3 GHz and 4 at 8.4 GHz for a total of 32 MHz in each 
frequency band. 

The data were correlated with the VLBA correlator at 
the Array Operations Center in Socorro, New Mexico. The 
correlated data were calibrated and corrected for time and 
frequency dependent phase variations usin g NRAO 's Astr o- 
nomical Imaging Processing System (AIPS. lGreisenl l|1998l) ). 
Initial amplitude calibration for each IF was accomplished 
using system temperature measurements and gain curves 
generated during observations. Fringe-fitting was done in 
AIPS using a solution interval equal to the scan duration 
and a point source model. 

Amplitude calibration was improved in a second stage 
by using observations of sources whose core flux density is 
known to be ^ 90 per cent of the total flux density. For a pre- 
cise definition of the core, refer to Section 3.1. To this end, a 
single amplitude gain correction factor was derived for each 
antenna based on fitting a Gaussian model to the core com- 
ponent. Gain correction factors were then calculated based 
on the difference between the observed and the model visi- 
bilities. Finally, the amplitude gain correction factors were 
applied to the target sources. This is a non-standard pro- 
cedure that improves the overall amplitude calibration. The 
accuracy of the amplitude calibration determined in this way 
is conservatively estimated to be within 20 per cent. 

Data imaging and modelfitting were done us- 
ing the Caltech Difference Mapp ing program 'Difmap' 
l|Shepherd. Pearson fc Taylor! 1 19951 ). after inspecting the 
data and editing out obvious bad points. The data were im- 
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Figure 2. Two examples each of class 'A', 'B', and 'C sources, respectively. Our current sample contains no examples of classes 'D' and 
'E'. 
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Tabl e 1. Optical and ra dio properties of the sample. The positions shown above are the most recent ICRF positions 
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aged using Difmap in automatic mode. Generally, this mode 
fails for about one third of the sources which have structure 
that is complex or too extended for the automatic script 
to handle. These have to be redone by hand in interactive 
mode. For the southern sources in this sample, almost all 
had to be imaged by hand due to the poor u-v coverage. 

Difmap combines the visibilities in all four 2.3 GHz re- 
ceiver IFs but does not correct for spectral index effects. 
It was assumed that the source structure variations across 
the IFs were negligible (6% variation in frequency at this 
band). Uniform weighting was used for the initial phase self- 
calibration before changing to natural weighting. For our ar- 
rays, uniform weighting gives more weighting to the longer 
baselines whilst natural weighting gives more weighting to 
the shorter baselines. Images of six of the 31 sources in our 
sample are shown in Figure [5] with the remainder available 
online. Images of additional epochs can be found in the form 
of contour plots at the RRFID and BVID websites. 

Generally, circular Gaussian models were used to fit the 
u-v data in order to parametrise the source morphology. 



Like imaging, modelfitting is an iterative process. Elliptical 
Gaussian components were used only to represent the core 
component or a very bright jet component if the residuals re- 
maining from a circular Gaussian model were too large and 
made it difficult to continue modelfitting using the residual 
map. The modelfits generally describe the visibility data well 
but these models may not be unique because of incomplete 
sampling in the u-v plane. In order to determine suitability 
of a source as a calibrator, we then determined the amount 
of source structure, as well as its variation with time, using 
several different methods based on both the CLEAN compo- 
nents from the Difmap imaging as well as on the parameters 
of the fitted Gaussians. 



3 ANALYSIS AND RESULTS 

Characterisation of the morphology of an AGN is very com- 
plex and no single metric is able to adequately define a good 
calibrator. Thus we have developed a number of complemen- 
tary approaches that are described below. A weighted com- 
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2 


30 


2 


05 


0.89 


6.58 


3.49 





53 



Table 2. Results of variability analysis in the flux density, core fraction and radial extent. 



bination of all of these metrics is eventually used to classify 
the suitability of each source to be a high resolution phase 
calibrator. These are discussed below. 

3.1 Core Flux Density 

A good calibrator should be relatively bright at the fre- 
quency of observation to be easily detectable. It should also 
be stable with minimal flux density variation over time. 
There is no precise definition of the AGN 'core' in the liter- 
ature. In general the bright, compact flat spectrum feature 
is referred to as the core. We have confirmed our identifica- 
tion of the cores for our entire sample by establishing that 
these structures have flat or inverted spectra. This was done 
by using the 8.4 GHz data that are observed simultaneously 
with the 2.3 GHz observations presented here. Here we de- 
fine a core flux density, S CO re for each source and for each 
epoch as the flux density of the Gaussian component fitted 
to the core. The latest core flux density as defined in this 
section is shown, for each source in Table [T] From the core 
flux density, we computed the mean core flux density (S) 
averaged over all epochs in which the source was observed. 
The extent to which the core flux density varies over time 
can be characterised by the core flux density variability in- 
dex (ffcore/S) where a core is the core flux density standard 
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Figure 3. Distribution of the core flux density variability index. 

deviation. A value of 0.0 indicates no variation over time 
(Table H column 5). 

Figure [3] shows that 90 per cent of the sources have 
a core flux density variability index below 0.4 and 23 per 
cent have a variability index below 0.2 (which is below our 
estimated calibration uncertainities) . These low variability 
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Figure 4. Distribution of the core fraction variability index. Note 
that all sources have an index ^ 0.4. 
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Figure 6. Distribution of the Weighted radial extent variability 
index. Seven sources with variability index greater than 1 are not 
shown in the histogram. 
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Figure 5. Distribution of the core fraction. This measure was 
used instead of the variability index as it shows slightly higher 
variation than the core fraction variability index. 



indices indicate that most core fluxes in this sample are sta- 
ble. 



3.2 Core Fraction 

Following lOiha et al.1 (|2004bh we define a core fraction as: 

C = S coreu2c / S total (1) 

where S CO re_cc is the sum of the flux densities of CLEAN 
components within one synthesised beam of the brightest 
pixel and Stotai is the sum of all the CLEAN component 
flux densities. 

It provides an indication of how point-like a source is 
and also provides a way to track source structure changes 
from epoch to epoch. The average core fraction, C, was com- 
puted for each source over all epochs. As with the core flux 
density variability index, we also computed the source core 
fraction variability index (ac /C) where ac is the core frac- 
tion standard deviation. A value of 0.0 for the source core 
fraction variability index indicates no variation over time. 
The results are shown in Table [2] column 8. 



The distribution of the core fraction variability index 
is shown in Figure [4] while the distribution of the core frac- 
tion is shown in Figure [5] The mean core fraction for all 
the sources is 83 per cent with a standard deviation of 12 
per cent. All sources have a mean core fraction variability 
index below 0.40. Twenty seven sources have a variability 
index between and 0.2 and the remaining 4 between 0.21 
and 0.40. In general, the southern sample sources are very 
compact with little variation. For this reason, we used the 
actual core fraction and not the variability index to classify 
the sources. 



3.3 Flux Weighted Radial Extent 

The previous two metrics are primarily measures of core- 
dominance. This and the following metric provide comple- 
mentary information by quantifying the radial extent of the 

sources in the Southern sample. 

The flux density weighted radial extent l|Oiha et al.l 

2004b) is defined as: 



R\ 



Et Sin 



(2) 



where Rw is in units of milliarcseconds and n is the radius 
at which the ith CLEAN component has flux density Si. 

The mean (Rw) an d standard deviation were also cal- 
culated as was the variability index (<jr w / Rw) where <tr w 
is the weighted radial extent standard deviation, see Tabled 
column 11. 

The weighted radial extent variability index is widely 
distributed as shown in Figure [6l with 10 per cent of the 
sources having an index between and 0.2, showing minimal 
variation in the weighted radial extent. 32 per cent of the 
sources have an index between 0.21 and 0.40 while 13 per 
cent have an index between 0.41 and 0.60. 16 per cent of the 
sources have an index between 0.61 and 0.80 while the rest, 
(29 per cent) have an index greater than 0.81. This shows a 
high degree of variability in the weighted radial extent. 
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Figure 7. Distribution of the Unweighted radial extent variabil- 
ity index. Six sources with variability index greater than 1 are 
not shown in the histogram. 

3.4 Unweighted Radial Extent 

This is the radial extent within which 95 per cent of the 
source flux density is contained. This provides a measure 
of how extended a source is. The mean (Ruw) and stan- 
dard deviation were calculated as was the variability index 
(<tr uw I Ruw) where ur uw is the unweighted radial extent 
standard deviation and Ruw is in units of milliarcseconds, 
(Table [H column 14) 

Like the weighted radial extent variability index (see 
above), the unweighted radial extent variability index is also 
widely distributed. Figure [7] shows that 6 per cent of the 
sources have an index below 0.2, 23 per cent have an in- 
dex between 0.21 and 0.40 while 29 per cent have an index 
between 0.41 and 0.60. About 10 per cent have an index be- 
tween 0.61 and 0.80 while 35 per cent have an index greater 
than 0.80, showing very high variability in the unweighted 
radial extent. As discussed below, this high variability (in 
both the weighted and unweighted radial extent) is largely 
a result of low brightness 'jet' features that are detected in 
some epochs but not others, usually due to differing sensi- 
tivities of the observing array rather than a real variation in 
the source. 



4 SOURCE CLASSIFICATION 

To capture the complexity of determining whether a radio 
source is a suitable phase calibrator for a radio interferome- 
ter, we have developed the set of complementary metrics de- 
scribed above. Here we discuss how we combine the informa- 
tion these metrics provide to arrive at a classification scheme 
for potential calibrators. We explain our reasons for weight- 
ing the different metrics as we have, and present our rec- 
ommendations based on this weighting scheme. Since such 
a scheme is necessarily somewhat subjective it should be 
considered as reasonable and useful rather than a definitive 
classification. 

The four metrics we used to classify the calibrator 
source quality were core fraction, core flux density variabil- 
ity index, weighted and unweighted radial extents. Each of 
these metrics is first assigned a score, with the scores as 



Core flux density 


Score 


Core fraction 


Score 


(%) 




variabilty Index 




> 1 


0.00 






0.91-1.00 


3.5 


91-100 


35.0 


0.81-0.90 


7.0 


81-90 


31.5 


0.71-0.80 


10.5 


71-80 


28.0 


0.61-0.70 


14.0 


61-70 


24.5 


0.51-0.60 


17.5 


51-60 


21.0 


0.41-0.50 


21.0 


41-50 


17.5 


0.31-0.40 


24.5 


31-40 


14.0 


0.21-0.30 


28.0 


21-30 


10.5 


0.11-0.20 


31.5 


11-20 


7.0 


0.01-0.10 


35 


1-10 


3.5 



Table 3. Source core fraction and core flux density variability 
index score distribution. 



Index 


Score 


> 1 


0.00 


0.91-1.00 


1.5 


0.81-0.90 


3.0 


0.71-0.80 


4.5 


0.61-0.70 


6.0 


0.51-0.60 


7.5 


0.41-0.50 


9.0 


0.31-0.40 


10.5 


0.21-0.30 


12.0 


0.11-0.20 


13.5 


0.01-0.10 


15 



Table 4. Distribution of scores for the weighted and unweighted 
radial extent variability index. 

shown in Tables [3] and [4] The scores for the individual met- 
rics were chosen so as to sum to 100 for a perfect calibrator 
source. We also chose to score the core fraction and core 
flux density variability so as to give them higher weight, 
each having a maximum score of 35, while the two radial 
extents have maximum scores of 15. We use a lower maxi- 
mum score for the radial extents as they are highly sensitive 
to small epoch-to-epoch variations in the signal-to-noise ra- 
tio of the images. Faint extended features may be detected 
only at some epochs depending on the signal-to-noise, re- 
sulting in a large variation in the radial-extent measures. 
Such low surface brightness features generally have negligi- 
ble impact on the usefulness of the source as a calibrator. 
In any case, for a lower resolution array like meerKAT these 
extended structures will be embedded in the main central 
component. 

Based on the overall score summed up from this weight- 
ing scheme each source falls into one of five classes 'A' 
through 'E'. A source falls into class 'A' if its overall score is 
between 80 and 100 per cent, class 'B' if the score is between 
60 and 80 per cent, class 'C if it scores between 40 and 60 
per cent, class 'D' for any score between 20 and 40 per cent 
and class 'E' for sources scoring between and 20 per cent. 
We propose the following classifications: 

A - Excellent calibrator (score between 80 and 100). 
B - Very good calibrator (score between 60 and 80 ). 
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Source Score Total Score Class 

Name Wi W 2 W 3 W 4 % 

A 
B 
B 
B 
A 
B 
B 
A 
B 
B 
B 
B 
B 
A 
A 
B 
A 
B 
B 
B 
B 
A 
B 
B 
A 
B 
B 
C 
C 
A 
C 



Table 5. Classification of the Sources. Wi - Core Flux Den- 
sity, W2 - Core fraction, W3 - Unweighted Radial Extent, W4 - 
Weighted Radial Extent 

C - Good calibrator (score between 40 and 60). 
D - Use with caution (score between 20 and 40). 
E - Unsuitable as calibrator (score below 20). 

In our sample (Table [S]), we have 9 class 'A' sources, 
19 class 'B' sources and 3 class 'C sources. There are no 
sources in class 'D' and class 'E'. Thus all the sources in 
our sample would be 'Good' calibrators and all but three 
sources are likely to be 'Very Good' calibrators. 



5 CONCLUSIONS 

We modelfitted up to 32 epochs of observations (average of 
20 epochs) for each of the 31 sources in our samp le which was 
select ed from the RRFID Kinematic Survey (|Piner et al.l 
120071 ) and determined their suitability as phase calibrators. 
While the kinematic survey looks more at proper motion in 
the sources at 8.4 GHz, this paper concentrates on the mor- 
phological properties of the sources at 2.3 GHz, a frequency 
more relevant to emerging southern hemisphere arrays like 
meerKAT and ASKAP. 

We have developed a method to classify radio sources 
according to their suitability as phase calibrators for radio 
interferometers. We first characterize a source by calculating 



several metrics which give measures of the degree to which 
the core dominates the source, and the degree of variability, 
both in flux density, and the degree to which the source 
is extended. These metrics are then combined to give the 
source a total score, which is used to assign the source to 
one of five classes of suitability as a calibrator. All 31 sources 
in our sample were classified as 'Good' calibrators with 28 
classified as 'Very good' or better. 
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